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Abstract 

We analyse a supersymmetric mechanical model derived from (1 + l)-dimensional field 
theory with Yukawa interaction, assuming that all physical variables take their values in 
a Grassmann algebra B. Utilizing the symmetries of the model we demonstrate how for a 
certain class of potentials the equations of motion can be solved completely for any B. In a 
second approach we suppose that the Grassmann algebra is finitely generated, decompose 
the dynamical variables into real components and devise a layer-by-layer strategy to solve 
the equations of motion for arbitrary potential. We examine the possible types of motion 
for both bosonic and fermionic quantities and show how symmetries relate the former to 
the latter in a geometrical way. In particular, we investigate oscillatory motion, applying 
results of Floquet theory, in order to elucidate the role that energy variations of the lower 
order quantities play in determining the quantities of higher order in B. 
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1 Introduction 



Classical supersymmetry sets out to extend the unified treatment of bosonic and fermionic 
quantities in the usual QFT framework to the classical level. Normally, in semiclassical treat- 
ments the fermionic variables are set to zero as soon as the supersymmetric theory has been 
constructed. The usual argument goes that since we cannot find classical fermions in nature, 
fermionic quantities should be omitted altogether at the classical level. 

However, this is far from necessary. In fact, a consistent approach to classical supersymmetry 
has long been available - for a review see e.g. the book by de Witt 0. Fermionic quantities are 
then treated as anticommuting variables taking values in a Grassmann algebra B. Grassmann- 
valued mechanics has been analysed in the works of Berezin and Marinov [|I| and Casalbuoni 
and later by Junker and Matthiesen p. A main difference to our work is that both [|l| 
and do not distinguish clearly between generators of the algebra and dynamical quantities 
and thus define the Grassmann algebra B rather implicitly. The fact that the bosonic variables 
take values in the even part of the same algebra B is not apparent in these works, although 
both recognize that the bosonic variables cannot be real functions anymore - without, however, 
elaborating on this fact. A central aim of this paper is therefore to make sense of the general 
Grassmann-valued equations of motion, including the fermionic ones, and to find ways to their 
solution, which is done in |I| and only in very special cases. Junker and Matthiesen, who 
investigate a similar mechanical model, achieve a more general solution than in and [Q], 
but again under the (implicit) assumption that the Grassmann algebra is spanned by only two 
generators identified with the fermionic dynamical variables. We can confirm most of their 
results (in different form, though, due to a different choice of variables) as special cases of our 
solutions. However, we disagree about some details, in particular, concerning the case of zero 
energy. 

The mechanical model that we study here is the supersymmetric motion of a particle in a 
one-dimensional potential, derived by dimensional reduction from the usual N = 2 supersym- 
metric (1 + l)-dimensional field theory with Yukawa interaction. A slightly different version of 
this model was investigated in 0, where a different concept of reality was used that led to a 
negative potential in the Lagrangian. The approach taken here stays closer to the usual case 
with the positive potential. 

An important result of [0 was that a complete solution for the particle motion could be 
found on the assumption that the underlying algebra B has only two generators. This led to 
relatively simple results, however is unnecessarily restrictive. 

Here we show first that for a large class of potentials the solution to the equations of motion 
can be found for any B and depends only on a small number of i3- valued constants of integration, 
one of which is a Grassmann energy E. 

To deal with essentially arbitrary potentials we adopt a second method which is closer to 
that of 0, although we need not restrict ourselves to two generators: Choosing the Grassmann 
algebra to be finitely generated, with n generators, we split all dynamical quantities and equa- 
tions into their real components, named according to the number of generators involved in the 
corresponding monomial. Then, beginning from the zeroth order equation, which can be seen 
as a form of Newton's equation, we subsequently work our way up to higher and higher orders, 
utilizing the solutions already found for the lower levels. This layer-by-layer strategy allows us 
to solve the equations of motion for any potential with reasonable mathematical properties. 

The existence of a complete solution to the coupled system of equations of motion looks 
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surprising in view of the increasingly large number of equations involved for large n. However, 
on second thoughts it is not so unexpected: Due to our first solution method we know that a 
full Grassmann solution can be found in many cases, the decomposition of which should give 
us exactly the component solutions obtained by the second method - which it does indeed as 
we shall demonstrate. 

A final word has to be said about the assumption of only a finite number of generators since 
it has been claimed that this must necessarily lead to contradictions: Emphasizing that our 
paper deals with the classical theory we do not find this to be true. 

We begin our analysis in section |^ by presenting the Lagrangian and the equations of 
motion that we will be concerned with in this paper. Essential for solving these equations are 
the symmetries and associated Noether charges of the Lagrangian which we therefore examine 
in section ^. For a certain class of potential functions, namely those for which a particular 
integral can be calculated analytically, we describe in section § how the equations of motion 
can be solved completely and illustrate this method for two exemplary potentials, the harmonic 
potential U{x) = kx and the hyperbolic potential U{x) = ctanhA;a;. 

Sections ^ and ^ are devoted to the description of our layer-by-layer method which is then 
explicitly carried out up to fourth order, and illustrated by the harmonic oscillator case in 
section |0| . 

Next, we investigate the symmetries in component form in section |^: While all component 
charges can be simply derived by decomposing the original charges, they also reflect a huge 
number of symmetries of the highest order component Lagrangian by which they can be found 
using Noether's procedure. In addition to the symmetries known from the original Lagrangian 
there appear extra classes of symmetries which cannot be related to the former. It is possible to 
trace them back to invariance under certain changes of base in the Grassmann algebra, viewed 
2"-dimensional vector space. 

One question that immediately arises from the decomposition of the physical quantities 
into components is addressed in section |^: How must we interpret the numerous new functions 
that arise from this procedure? Utilizing the new symmetries found in section |^ we offer a 
geometric interpretation for the lowest order bosonic and fermionic functions and describe the 
three basic types of motion possible for them. The higher order quantities are interpreted here 
from a different point of view: We see them as variations of the lower order quantities with 
respect to the integration constants involved in these functions, namely an initial time Iq and 
the energy Eq. 

In section ^ we apply our results to study general oscillatory motion, by which we mean that 
the lowest order bosonic function is periodic. The characteristic appearance of linear-periodic 
terms is explained from both a physical and mathematical perspective employing results of 
Floquet theory. Again, we choose a particular potential to illustrate this in section |9.1| , which 
allows us also to show that the two solution methods presented in this paper coincide. 

A common restriction on both methods is that the energy Eo of the system must be positive, 
therefore the zero energy case has to be discussed separately. We do this in section p!o| . 

Some ideas for further generalization and analysis conclude this paper. 
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2 Super symmetric Mechanics 



We start our discussion with the standard Lagrangian density for (1 + l)-dimensional super- 
symmetric field theory with Yukawa interaction: 

I I i i dU 

C = -d+cj)d^(t) - -U{(f)f + d-ip+ + -tp- d+tp- + i—ip+i)-, 

containing a real bosonic scalar field 0, a real two-component fermionic spinor field ijj and a 
potential function U{(f))] d± are the light cone derivatives dt ± dr,.. 

We assume that (j) and the two components of ip take their values in the real even and odd 
part, respectively, of an arbitrary Grassmann algebra B, and that the potential function U{(p) 
can be expanded into a power series in with real coefficients. Complex conjugation is defined 
such that {zi Z2)* = 22-2* in accordance with the conventions in so that the Lagrangian 
density is a real function despite the z-factors that occur in front of the fermion terms. As 
we are dealing with the classical case we assume further that all fields and their derivatives 
commute or anticommute, depending in the usual way on the bosonic or fermionic nature of 
the fields. 

In this paper we shall be interested only in spatially independent fields, i.e. dx4> = dxi}+ = 
dxip- = 0. As the fields are then functions of time only this leads us directly from field theory 
to mechanics. It is therefore sensible to think of the bosonic field as describing the one- 
dimensional motion of a particle in the potential [/^. To support this notion we will change the 
variable to x for the rest of this paper. We then obtain the following Lagrange function: 

L = \i^- \u{xf + + + I U'{x)^+ij. . (1) 

where the dot denotes a time derivative and U' means the derivative of U with respect to x. 
Performing a formal variation of this Lagrangian with respect to the variables x, -0+ and ifj- 
and neglecting total time derivatives we derive the equations of motion for the system: 

X = -U{x)U'{x) +tU"{x)ilj+-^_ (2) 
^+ = -U'{x)^p_ (3) 
ip. = U'{x)il)+. (4) 

There is a slight ambiguity in these equations if the Grassmann algebra is finitely generated 
and has an odd number of generators. Then the two equations for the fermion variables will 
be determined only up to an arbitrary function of highest order in the Grassmann algebra. We 
can think of this as a gauge degree of freedom and will come back to this point later. 

In order to solve the equations it is advisable to understand the symmetries of the system 
first. This we will do in the next section. 



3 The symmetries of the model 

The first thing to notice is that the Lagrangian has no explicit time dependence. We therefore 
have invariance under time translation, leading to a conserved Hamiltonian as the corresponding 
Noether charge. This we calculate to be 

H = \i^ + \u{xf - iU'{x)^+ij_; (5) 
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its conserved value we call the energy, denoted by E. Note, however, that this Hamiltonian 
is an even Grassmann-valued function and that therefore the conserved energy E is an even 
element of B. 

The operator ^ generates time translations, and acts on the dynamical variables in the 
obvious way: 

Ax = r]x, A-01 = r]iiji, = (6) 
where 77 is an infinitesimal even Grassmann parameter. 

In addition to time translation invariance there are two further independent symmetries of 
the Lagrangian, relating fermions and bosons. These can be written in infinitesimal form as: 

Sx = ^e'0+, Sip+ = —ex, Si/j^ = —eU{x), , , 

6x = ieip-, 6ip+ = eU{x), 5ip- = —ex, 

where e is an arbitrary infinitesimal odd Grassmann parameter. These transformations lead 
only to a change in the Lagrangian by a total time derivative. We can therefore apply Noether's 
procedure and find the following charges: 

Q = xilj+ + U{x)ij_, (8) 
Q = xilj.-U{x)i{j+, (9) 

the conservation of which can be easily shown using the equations of motion. Note that the 
charges are odd elements of the Grassmann algebra. 

As for time translation invariance we now define two operators Q and Q generating the 
two symmetry transformations. From we can read off the action of these operators on the 
dynamical variables: 

Qx = ip^, Qip+ = ix, Q V- = iU{x) 
Q X = ijj-, Q 'ip+ = —iU{x), Qijj_ = ix. 

Using the action of the operators on x, -0+ and we find that Q, Q and ^ form a closed 
algebra with the relations 

as long as the equations of motion are satisfied. Notice that ^ commutes with everything. 

This is the usual N = 2 supersymmetry algebra in space dimensions, and with (^) we have 
an on-shell representation of this algebra. Thus we will from now on speak of supersymmetry 
transformations and call the associated charges supercharges. 

Beside time-translation invariance and the two supersymmetries there exists still a further 
invariance, this time for the fermionic functions only. In infinitesimal form it is given by the 
transformation 

= r7^_, Aip_. = -r]ip+, (11) 

where r] again denotes an infinitesimal even Grassmann parameter. We can think of this as 
an internal rotation of the fermionic variables. This invariance leads to a further conserved 
Noether charge: 

R = iip+ilj_. (12) 
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Comparison with (|l]) shows that the fermionic functions enter the interaction term in the 
Lagrangian only via R. This leads directly to the fact that the x-motion depends on the 
fermionic functions only through this one constant. 

In addition, the Hamiltonian reveals that the only fermion contribution to the total 
energy E is through R. This is an important simplification, and it will allow us to solve the 
equations of motion completely for a number of potential functions U without restriction on 
the nature of the Grassmann algebra B. 



4 General solutions via Grassmann integrals 

We have already mentioned that the x-equation of motion (|]) nearly decouples from the other 
two equations because the fermion functions in the coupling term i U'{x) form a Grass- 

mann constant. It is therefore sensible to begin with this equation. 
From the conserved Hamiltonian we know that 

= 2E -U{xf + 2RU'{x). (13) 

For the next step we have to use that every Grassmann number z can be split into two parts, 
its 'body' and its 'soul': z = Zb + Zg. The body is just the real number content of the soul 
is the remaining linear combination of products of (odd) Grassmann generators and will be 
nilpotent if the Grassmann algebra is finitely generated. Note, however, that we do not need 
such a restriction for body and soul to be well-defined. A square root for a Grassmann quantity 
can be defined by its power series as long as the body is positive. Since R is the product of two 
odd Grassmann terms, its body and therefore the body of the whole third term on the right 
hand side of ([T^) is zero, leaving the restriction 2Ei, — U^{x) > if x is to be well-defined. This 
just means that the kinetic energy of the classical particle moving in the potential U"^ has to 
be positive. 

The resulting first order differential equation for x is 

dx 



±^2E -U{xY + 2RU'{x). 

Provided B is finitely generated, we may regard x as lying in the vector space B and apply the 
standard theory of systems of ODE's to show that this equation has a unique solution for any 
given initial data a;(to)- 

The equation can formally be solved by separating variables: 

t-to = ± . (14) 

■J^ito) ^2E - U{xy + 2RU'{x) 

Note that while the left-hand side is just a real expression, the right-hand side is a Grassmann 
integral of a Grassmann integrand. Such an integral is defined as a line integral in the Grass- 
mann algebra, thought of as a finite or infinite-dimensional vector space spanned by products 
of generators. As is shown in |^ such an integral is independent of the actual path relating 
start- and endpoint of the integral since the integrand is Grassmann-even. 

Now the integrand comes from an ordinary real function extended into the full Grassmann 
algebra using its power series. We can therefore sensibly ask whether there is an indefinite 
integral F{x) to this real function. If so, we can extend this function back into the Grassmann 
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algebra thus gaining an indefinite integral for the Grassmann integrand. Because the integrand 
is even it then follows that the integral is given by the difference of F evaluated at the start- 
and endpoints of the path: 

t-t^ = ±{F{x{t))-F{x{to))) 

If the function F has an inverse, we finally get a; as a function of t and therefore the solution 
of (^. We will illustrate this method below for two potential functions U . 

The solution to the fermion equations @) and @ is now easy. From the solution for x{t) 
we can immediately calculate x{t) and U{x{t)). Using the explicit formulae (||) and for the 
two conserved supercharges Q and Q we can evaluate the linear combinations 

Qx-QU{x) = 2Hij+, 
QU{x) + Qx = 2Hij_, 

where we have used that -R'0+ = Rip- =0 which follows from (p!2D . Since the Hamiltonian H 
has the constant value E we deduce that 

= (15) 

4,_ = om^Qt, (16) 

This effectively solves and as the values of Q and Q are determined by the data at the 
initial time to. There is one subtle point here, however, since the division by 2E is only possible 
as long as E^, the body of E, which can be interpreted as the classical energy of the particle, 
is non-zero. Indeed, when we analyse the solution for finitely generated Grassmann algebras 
later, we will find again that the (i?;, = 0)-case is special. 

A second point which has to be mentioned here is that we have treated R as an independent 
parameter. However, R is determined by the two fermionic quantities ■0+ and ip- and, in effect, 
by the two supercharges Q and Q: 

R = « 77, 

2E' 



which can be verified using the definition of R in (|T2|). This does not render our solution invalid 
but it shows that R is not a parameter that can be chosen independently. But let us now turn 
to some examples: 



4.1 The harmonic potential U{x) = kx 

We start with one of the easiest problems, the harmonic oscillator with potential U{x) = kx, 
k being real. Here the integrand in (|I^) is given by: 

1 

^{2E + 2kR) - {kx f 

where E is the constant Grassmann energy. Note that due to the special nature of the potential 
we can combine the U' {x)-ieim. 2kR with the energy E into one overall constant. If we treat 
the integrand as an ordinary real function it has an indefinite integral: 

/ / o ^ ,9 o dx=\ arcsin \-x \ , 
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so that we can write 

x{t) 

I fi: 

t-to = ± 



1 • ^ 

— arcsmi , = x) 

k W'2E + 2kR ' 



This formula can be inverted to yield the solution 



xit) = x(to) COS k(t - to) ± sin k(t - to), (17) 

k 



where we have denoted the constant + 2kR — k'^x{to)'^ by v{to)- This looks formally like 
the usual harmonic oscillator solution, and it is periodic with period but all terms are 
Grassmann- valued and not just real functions. 

For the two fermion terms we need only calculate x and U (x) and insert them into equations 



(p!5|) and (p!6D, which leaves us with: 



-0+ = p+ cos k{t — to) ~ p- sin k{t — to) , (18) 
i/j^ = cos k{t — to) + p+ sin k(t — to), (19) 

where 

_ ±Qvito) - Qkxjto) _ ±Qvito) + Qkxjtp) 
~ 2E ' ^' ~ 2E • 

With (plTl), (pT^) and (|19D we have found the complete solution to the equations of motion for 

Eh > 0, independent of the nature of the Grassmann algebra B. 

4.2 The hyperbolic potential U{x) = ctanh/cx 

As a second example we choose the hyperbolic potential U{x) = ctanhfcx, with c and k both 
real. The integrand in ( p!^ is 



{2E + 2Rck) - (c2 + 2Rck) tanh^ kx 



and - viewed as a real function - has an indefinite integral, the form of which depends on the 
constants involved: 



-ji — arcsin ( — 1 sinh x] : - > 1 
I =dx = - <^ ^arcsinh (Jl - ^ sinhx) ; - <l (20) 



a-6tanhfcx " iginhx 

Because the Grassmann-valued function is completely determined by the corresponding real 
function we have to look at the body of 

h c2 + 2ckR 



a 2E + 2ckR 

to decide which integral to use. Since R as product of two odd Grassmann constants has no 
body, the crucial term is 
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All three integrals (pO|) can be inverted to give x as a function of t: 

^ r arcsinh(y2^ysin(^,(t_to) + ■,cl>2E, 
^W = ^ ^rcsmhl^^^smh{ujjit-to) + Kn)) ■,cl<2E, (21) 
[ arcsinh (a;///(t-to) + /tm) ; = 2Eb 



where 



= k^/c^-2E , Ki = arcsin {\f- 2E+2ckR sinh(A;x(to)) 
ujn = ky/2E-c^ , Kjj = arsinh (y^^^^^sinh(A;x(to)) 
i^iii = k\/2E+2ckR , Km = sinhA;x(to)- 

From equations (plSf ) and (^6|) we can calculate the fermion solutions: 



^+ I ^ 1 V2ii;f2cfcfl I iQ-QVai^^ta'^h?/^^ ■c'^<2E (22) 

^- J \ ±g + gy^tanhy,, ' " 

i/j- j 2EVi^ \ ±gV2^+2cA;/? + Qcym ' ' 

with yi = ±(jji {t — to) + Kj, etc. Again, we have found the complete solution for the system 
without using any information about the Grassmann algebra B itself. We will return to this 
example later, after we have investigated a special class of Grassmann algebras. 

5 Finitely generated Grassmann algebras 

We will now make a choice for the underlying Grassmann algebra B of the system, namely that 
it is generated by a finite number of elements = 1,. . .,n with the property 

Every element z oi B can be written in the form 

n 

z = ^ ■Zji...ifc^n ■ ■ ■ 
k=0 

with complex coefficients 2^11. .ij.. zq is the body of z. The requirement for z to be real (in 
the sense explained earlier) fixes each coefficient to be either purely real or purely imaginary, 
depending on the number of generators involved. For practical reasons we will usually restrict 
ourselves to the case of four generators, although the general aspects of our discussion should 
be true for an arbitrary number n. 

In the case of four generators it is useful to define the following real monomials: 

6j = ^ijk = i^i^j^k, '■ ^1234 = — ^lS^3^4- 

Note that due to antisymmetry there are only six linearly independent monomials involving 
two generators, four monomials involving three generators, and just one monomial of highest 
order, i.e. involving all four generators. 
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We can decompose the dynamical quantities x, i/j^ and iIj^ as follows: 



i)-{t) = Pi{t)^i + Pijk{t)iijk, 

where it is implied from here on that the summation over indices is with i < j < k. In total 
there are 24 independent real functions which we call of zeroth, first, second, third or fourth 
order according to the number of generators involved in the corresponding monomials. 

Every analytic function involving x, ip^ or can be decomposed as well using its Taylor 
expansion. Applied to the potential function U this yields: 

U{x) = U{xo) + U'{xo)XijCij + (^f/'(xo)xi234 + ^U" {xo)x[i2X34.]j ^234, 

where the brackets denote antisymmetrization. For example 

a;[12a;34] = X12X34 — X13X24 + 3^142^23 + a;232;i4 — 2;24a;i3 + Xs4Xi2- 

Note that U{xo) and its derivatives are ordinary real functions. Inserting these results into (|l|) 
yields eight component Lagrangians of order zero, two and four, respectively: 

Lo = 1^0' -lu' (23) 

Lij = xoXij - UU'xij + XiXi + piPi + U'{XiPj - piXj) (24) 
1 1 

^1234 = -^[122:34] + a;oil234 - - (UU')' X [12X34] - UU' Xi234 (25) 
+ A[i23A4]+/9[i23P4] +U' (A[i23P4] - P[123 A4] ) + U" X[i2 ( A3/94 - A4P3)] 

where the argument of U, U' and U" is always xo{t), here and below. 

From the highest order Lagrangian L1234, which is a functional of 24 generally different 
component functions, we get the following set of Euler-Lagrange equations: 

Xq = -UU' (26) 
x,j = -{UUjx,, + U'Wp, - pAj) (27) 

X1234 = -1{UU')" X[i2X34] - {UU'Y X1234 

+ U" (A[123P4] - P[123A4]) + U"'X[i2 (A3P4 " A4P3)] 



A. \ ^ / -U'pi 
U'Xi 



pi 



Xijk 1 _ I —U'pijk — U"x[ijPk] 
Pijk J 1 U'Xijk + U"x[ijXk] 



(28) 
(29) 
(30) 



The same equations can also be obtained by splitting the original equations of motion (^-(^) 
into their components. 

It is remarkable that all equations can be derived from just the one Lagrangian. The 
other seven component Lagrangians are completely redundant and yield no new equations: 
From (p3|) we can only obtain equation (p6D whereas from the Lagrangians of type (|2^) we can 
derive equations (|2|), (^) and (H). 
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Now we can also see where the ambiguity mentioned earher comes from when there is an odd 
number n of generators: Because the Lagrangian is an even functional, none of its components 
can contain any function of highest order n, so these functions cannot be governed by any 
equations of motion. E.g. for three generators the functions A123 and P123 will stay completely 
undetermined. Formally, we can derive evolution equations similar to ( |30D by decomposing (|^) 
and @ into components but they are meaningless since an arbitrary function can be added 
to both sides. Thus we can regard these highest order functions as non-physical and treat 
them as trivially separated gauge degrees of freedom. However, if there is an even number 
of generators this problem cannot occur, since there will always be one Lagrangian of highest 
order containing component functions of all orders which are thus determined by equations of 
motion. 



6 Layer- by- layer solutions 

We now proceed to derive the solutions to the equations of motion (|26|)-(|30D for an arbitrary, 
sufficiently different iable potential function U adopting a layer-by- layer strategy. That means 
that we will start with the lowest order or 'body' equation (|^) and then use its solution to work 
our way up to the higher order and more complex equations, including also the fermionic ones. 
We do this here up to fourth order, although there is no obstruction in principle to continue the 
strategy for an algebra with more than four generators. In fact, in every layer the equations 
and their solutions are the same irrespective of how many generators there are; there will only 
be a different number of them. 

Before we start with the bosonic equations, it is useful to look at the decomposition of the 
Hamiltonian H = Hq + Hij^ij + //1234C1234 and of the Noether charge R = Rij^ij + -R1234G234 
(note that R does not have a body): 

^0 = lil + lu'ixo) (31) 
Hij = XoXij + UU'xij — U' Rij (32) 

-f^l234 = ]jXyi2X^A\ + XqXi2M + ^ iPU')' X [12X^4:^ + UU' Xi2Zi (33) 
-[/'i?1234 - U"x[i2Rzi] 

Rij = XiPj - piXj (34) 

-R1234 = -^[123^4] "~ P[123'^4]- (35) 

Since these are components of Grassmann conserved quantities, they are conserved too, as can 
be checked easily using the equations of motion. All components can be derived as Noether 
charges from the highest order Lagrangian I/1234 and we will do this in section |^ but for now we 
only need that they are constant in time. The respective values of the Hamiltonian functions 
Ho, Hij and i^fi234 are denoted in the following by Eq, Eij and i?i234- 



We notice first that ( pGl) is just the standard Newtonian equation of motion for a particle 
moving in a potential ^U^- This is the bottom layer of the system. 

The solution to ( p6|) using the constancy of the Hamiltonian Hq is well-known: 

t — to = ± / dx'n. (36) 

^^o(to) ^2Eo - U^ix'o) 
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The sign has to be chosen carefully to comply with the direction of motion of the particle; 
if the particle motion changes direction, the integral will be multi-valued and has to be glued 
together from pieces with unique sign to ensure that the overall result for t{xo) is monotonically 
growing. The implicit function theorem then allows us to locally invert t{xQ) and obtain the 
required Xo{t). Since (^) is a second order equation there are two constants of integration, 
(^o) (or equivalently to) and the energy Eq. 

Whereas (p6D is a non-linear equation, (^) is an inhomogeneous linear equation (which 
can be simplified using (|5^). It is now convenient to use the Hamiltonian (^): Solving the 
equation Hij = E.ij for x-ij we can reduce the problem to a first order differential equation. This 
can be solved by standard methods to yield: 

CijXQ ~t~ Ej^jXQ I —^dt -\- RijXQ I -r2"t/t . 

Jto Xq Jto Xq 

where Cij is an arbitrary integration constant related to the initial value of Xij. The second 
integration constant of (pTj) is the energy variable Eij. This result was also found in with 
Rij set to —1. However, using equation (^) and Hq = Eq, the last term on the right hand side 
can be rewritten: 



Xij 



R,XQ -,dt = R,XQ -j^^dt = 1 -—^dt = ^U. (37) 

so that we end up with the simpler expression 

/■* 1 Ri 

Xij = CijXo + EijXo / ^dt' + ^r^U. (38) 
Jto xt^ 2Eq 

All three functions which occur in ( p8|) can be calculated directly from Xo(t). 

The same procedure can be applied to solve (pSD which simplifies when we use (^4]) and (|35|) . 
We rearrange the equation -^1234 = -E1234 to isolate i;i234 and then solve this first order equation. 
The result can be written as: 

2:1234 = Ci234io + -^1234*0 / -T^dt' + (39) 

Jto Xq ZhlQ 

+ Uq f "'"""J + Uq f . 

2 Jto Xq 2 J to Xq 

Again, there are two additional integration constants, C1234 and E1234. Note that the solution 
depends only on functions that we already know from the lower layers, namely Xq and Xij. In 
fact, the solution can be expressed as a function of Xq and its derivatives only by inserting (^) 
into (H): 



3^1234 



Cl234io + ^1234^0 / -Tndt' + ^^JT^U (40) 
Jto X^ 2Eq 

+ ^C[i2C34]Xo + C[12£^34] ( ^0 / -Todt' + ) 
2 ' ' ' \ Jto Xi Xq ) 



+ 2 ^[12-^34] 



2Eq 
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This complicated-looking expansion is particularly useful when we come to the interpretation 
of the solution in section R. 



There remain the fermion equations ( P^D and (^Up to be solved. Both are systems of two 
homogeneous first order equations and therefore have two independent solutions. Just as we 
used the decomposition of H to solve the bosonic equations, so it is now appropriate to de- 
compose the conserved supercharges Q = Qi^i + Qijk^ijk and Q = Qj^i + Qijkiijk to solve the 
fermionic equations. We find: 

Qi = XiXo + piU, 
Qi = Piio ~ XiU, 

Qijk = XijkXo + pijkU + X[iXjk] + P[iU Xjk], 

Qijk = Pijkio — XijkU + P[iXjk] — X[iU'Xjk]. 

Using suitable linear combinations of Xq and U with coefficients Qi and Qi and applying Hq = Eq 
we find that 

Xi = liXo - TiU, Pi = TiXo + kU, (42) 

where /j = and rj = hence the two integration constants of the solution are basi- 
cally given by the first order supercharges. A similar result was obtained in [^j, although the 
restrictions placed there on the shape of the potential U are unnecessary. 

Note, however, that this solution does not work if the particle energy Eo equals zero. This 
has to do with the fact that both xq and U must vanish in this case. We will return to this in 



section 10 



The same procedure can be used to find the solution to equation (^) and we get the result: 

A.,. = - nj,U - ^ + "'^""^ -^V'i„x,.,p,, ^^^^ 
P,. = r„,i, + l,,U - ^ - "'"'''^ - ^'^""W-^'l , (44) 

where we have used the equations Hq = Eq and Hij = Eij and the identities R[ijXk] =R[ijPk] = 0. 
The new integration constants Ujk and Vijk are related to the third order supercharges Qijk and 
Qijk in the same way as k and are to Qi and Qi. As was the case for the second and fourth 
order bosonic solutions we find that Xijk and pijk depend only on functions we know already 
from the lower order layers, namely, Xq, Xi, pi and Xij, and again we can simplify the solution 
even further, inserting (|38D and (|4^ ) into equations (|43|) and (^4]): 



/ ft 1 ^ 1 Xo\ 
Xijk = hjkXo + C[ijlk]XQ + E[ijlk] Xq / —dt' + — (45) 

\ Jto Xq Xq hin I 



I to Xq Eq ^ 



-TijkU - C[ijrk]U'xQ - EyijTk] {u'xq 
Pijk = rijkXQ + C[ijrk]XQ + E[ijrk] (^Xq ^dt' + ^ - (46) 



+kjkU + C[ijlk]U'xo + E[ijlk] i U'xq 



t y)- 

to Xq rjQ J 

We have now derived all solutions to equations (|26|)-(|30D, i.e. we have explicitly solved the 
equations of motion up to fourth order for an arbitrary potential function U . They are all 
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functions of only four quantities and their time derivatives: 



/■* 1 /■* 1 

Xq, U, / —dt' and / —dt'. 

J to Xq Jta Xn 



* 1 
to Xq 



To illustrate the method we choose the harmonic oscillator as our first example. This has the 
advantage that we can compare the solution with our earlier result (|1^-(|19]). 



directly, 
(47) 



6.1 The harmonic potential U{x) = kx 

The solution for xq could be found using ( P^ ) but here it is easy to solve equation 
with the familiar result: 



XQ{t) = xo(to) cosA;(t-to) + ^ sin A;(t-to); Vq = ±^2EQ-k'^xl{to). 
The functions xq and U can be calculated easily and we find: 



m = 
Mt) = 

where Li and Ri are the constants 

Ri 



Li cos A;(t — to) — Ri sin A;(t — tg), 
Li sin k(t — to) + Ri cos k{t — to), 



Qi 



Qi 



2K 



2^0 



-fcxo(to 
) 



(48) 
(49) 



(50) 
(51) 



Next we have to calculate x^ and thus we need 



* 1 

to 3^0 



dt' 



xo{t) 

XQ 



(to) 2i?o-Fx'o2^'^o 



1 (xo{t) XoitoY 



2Eo \xo{t) 



(52) 



On using (| 

Xij (t) 



we find that 
Ri 



2Ei 







kxo{tQ)+CijVQ ) cosk{t — to) + 



Ri E. 
vo-Cijkxo{to)-{ 



2E, 







kv. 



sinfc(t-to). (53) 



For Xijk and pijk we need not calculate any new terms except Xq. After collecting various 
constants into Lij^ and Rijk we end up with 

Xijk = Lijk cos k{t-to) - Rijk sin k(t-to) 
Pijk = Lijk sin k{t -to) + Rijk cos k{t -to). 

Finally we have to calculate 0:1234. Therefore we need the integral 



* 1 

to Xq 



dt' 



1 Xq ^ 1 Xo 



GEoxl GE^xoUoito) 



Xo{t) 



There seems to be a singularity when xq approaches zero but fortunately this term has to be 
multiplied by xo and combined with two other terms which share the same coefficient ii^[i2-E34]. 



13 



One finds that all terms involving xq in the denominator cancel each other, so that the overall 
result is just a linear combination of xq and xo, i.e. of cos- and sin-terms. This holds also for 
all other contributions to 0:1234, so that we can write: 

3:^1234 = Ci cos kit — to) + C2 sin k{t — tQ), 

where Ci and C2 are rather complicated functions of -E'1234, C1234 and all lower order integration 
constants. 

The most interesting observation is that all bosonic and fermionic functions are linear com- 
binations of cos- and sin-terms only. As we will see later, this similarity between the solutions 
of different orders is the exception rather than the rule and a special feature of the harmonic os- 
cillator. We also find that the fermionic functions oscillate with the same period as the bosonic 
functions and differ only in amplitude and phase. Once the motion of the lowest order bosonic 
function xq is fixed, the principal motion for Aj, pi is almost completely determined. This is a 
general feature of our system, and it can be fully understood by investigating the symmetries 
of our model. 

Finally, it is easy to verify that the complete component solution is compatible with the 



general solution found in section 4.1. 



7 Symmetries in component form 

We have already dealt with the decomposition of the Hamiltonian H, the constant R and the 
two supercharges Q and Q. Now we want to explain the origin of the component charges 
from the underlying component symmetries. The Lagrangian L1234 (^) is all we need to look 
at, since the non-trivial symmetries of all lower order Lagrangians form subgroups of the full 
symmetry group of -^1234- One would not guess the extraordinary number of symmetries of 
L1234 if one did not know their supersymmetric origin. 

To start with, is invariant under time translations, and this leads to the conservation 
of the highest order Hamiltonian ifi234. This result can be generalized by looking at the full 
Grassmann transformation (^. Decomposion into components gives us back the time transla- 
tion symmetry mentioned above when rj = r]o is real; the choices r] = rjij^ij (no summation) or 
V = ^1234 ^1234 lead to seven extra charges which correspond to the Hamiltonians and Hq, 
respectively. 

Apart from the Hamiltonians, sixteen Noether charges Qi, Qi, Qijk and Qijk derive from the 
supercharges Q and Q and belong to a set of symmetry transformations which can be obtained 
from by the choices e = ej^j and e = eijk^ijk (no summation both times), respectively. As 
an example, Qi and Qi are the charges that belong to: 

5iXi2M = , SjXjki = — ejio ? SjPjki = —^iU, 

'5ia;i234 = ^ipi 5 ^i^jki = ^iU , SiPjki = — Cjio) 

where {ijkl} is a cyclic permutation of {1234} and no summation over i is implied. The 
term supersymmetry transformations should be avoided here, though, since the dynamical 
quantities in component form are just real functions - therefore they can be termed bosonic or 
fermionic only by convention. Notwithstanding this fact the component transformations and 
the corresponding charges reflect the underlying supersymmetry. 
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Further component Noether charges derived from the Lagrangian (^) are Rtj and i?i234- 
The transformations can be read off from (^TJ) choosing r] = rjij ^ij (no summation) and 
7] = Tjo (real), respectively. In the latter case we get, for example 



(^1234^1 — VO Pi ; 5l234Pi — —VO^i, 
<^1234-^ijfc = VO Pijk j ^USAPijk = ~V0 \jk- 



(54) 



So far all transformations have just resembled those found in section |^. There are however 
further sets of symmetry transformations which only occur for the component Lagrangian L1234. 
The easiest of these is given by 

SiXjki = -Vi^i: ^iPjki = -ViPi: (55) 
where again {ijkl} is a cyclic permutation of {1234}. This leads to four charges 

S^ = l{Xl + pl)■ (56) 

The next easiest is 

SijXikl = TlijK-, SijXjkl = ~VijXji ^ijPikl = VijPiy ^ijPjkl = ~1ijPj (57) 

and leads to six independent charges 

Sij = XiXj + pipj. (58) 

We shall exploit these charges later but for now we are more interested in where they come 
from. It turns out that these and two other groups of transformations not mentioned yet 
reflect invariance of the original Lagrangian under change of basis in the Grassmann algebra. 
A Grassmann algebra with n generators can be viewed as a 2'^-dimensional vector space in 
the obvious way, but the choice of the generators and their products as basis vectors is 
somewhat arbitrary. A change of basis, however, has to be compatible with the multiplicative 
structure of the algebra. This means firstly that we only have to consider transformations of 
the n generators and secondly that we can change an (odd) generator only by another odd 
element of the algebra. The final constraint is that the highest order monomial, e.g. ^1234 in 
the case of four generators, remains unchanged by the transformation so that the highest order 
Lagrangian, here L1234, remains invariant. 

In the four generator case, this leaves only a relatively small group of acceptable linear 
transformations. There are eight independent odd basis elements, namely and ^ijk- Cor- 
respondingly, there are the following independent infinitesimal transformations that fulfil our 
conditions: 

1- 6^6 + V^j {i 7^ j) 

2. ^— > — Tj^jki {{ijkl} an even permutation of {1234}) 

3. ii + Tj^i, ^— — r]^j {i 7^ j): There are just three independent scaling transforma- 
tions of this type. 

4- ^ ii + vCiki, ^ — V^jki {{ijkl} an even permutation of {1234}): There are six 
independent Grassmann scaling transformations of this type. 
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In order for the Grassmann quantities x, i/j^ and '0_ to be invariant, corresponding to these 
transformations of the basis vectors of B there have to be transformations of the real components. 
It is these component transformations that we have found in (|55|) and (pTl), which belong to 
the basis transformations]^, and^., respectively. 

The component transformations belonging to |I[ and ^. are also symmetries of the La- 
grangian, but since they are more complicated to write down and we will not use them later 
we refrain from giving them here. It should be mentioned though that they generate a full SL^ 
symmetry group. 

8 Bosonic and fermionic motion — general properties 

We finally come to the physical interpretation of the results we have so far obtained concerning 
the component dynamical variables. Starting with the lowest order bosonic function xq we 
first recall that this always describes the motion of a particle in one dimension in the potential 
|t/^(xo). There are three possible types of motion in one dimension: 

1. Movement with no turning point: The particle velocity is always positive (or negative), the 
motion can be bounded, approaching finite values of xq as t goes to ±00, or unbounded. 
Example: Flat potential function U (x) = c. 

2. Movement with one turning point: The particle velocity changes sign once, when xq 
reaches a maximal (minimal) value. Again, the motion can be bounded or unbounded. 
Example: Reciprocal potential U{x) = ^. 

3. Oscillatory motion: The particle velocity changes sign infinitely often, thus the motion is 
always restricted to a finite interval. Example: Harmonic oscillator U{x) = kx. 

There is one conserved quantity, namely the lowest order Hamiltonian Hq, given in (0). We 
can interpret this as half the squared length of the two-dimensional bosonic vector 

( *° ) ■ (59) 

The conservation of Ho means that the motion of this vector is restricted to a circle with 
squared radius 2Eq. 

Corresponding to the three types of motion for xq there are three types of motion for this 
vector: 

1. Movement with no turning point: The vector moves on the right (left) semicircle in the 
xo-U-pla.ne. The number of direction changes depends on the number of potential extrema 
of U. 

2. Movement with one turning point: Motion starts on the right (left) semicircle of the xq- 
f/-plane, then changes when the turning point is reached to the other semicircle where it 
then mirrors the previous motion. 

3. Oscillatory motion: There are two subcases, depending on whether the number of sign 
changes of U (x) between the two extremal points Xmin and x^ax is odd or even. 
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(a) When there is an odd number the vector will continually wind around the circle; 
direction changes depending on the number of extrema of U can be superimposed. 
This is the only case where the motion covers the complete circle. The most impor- 
tant example is the harmonic oscillator where we can see directly the circular motion 
of the bosonic vector. 

(b) When the number of sign changes of U is even, the bosonic vector will oscillate on 
the circle between two points lying symmetrically on either side of the [/-axis. An 
example of this case is the potential function U (x) = \J {kxY + c with c > 0, which 
leads to a harmonic oscillator potential shifted upwards by c. 




The bosonic vector is important because it largely determines the behaviour of the n first 
order fermionic functions. To see this, we look at the n two-dimensional vectors 

: = l,...,n. (60) 

The charges Si introduced in (|56|) guarantee that the lengths of these vectors are conserved, 
thus their motions are also restricted to circles. The squared radii are given by 25'^, providing us 
with a nice geometrical interpretation of these constants. Furthermore, the charges S'jj, given 
by (0), can be interpreted as scalar products between the zth and jth vectors, thus effectively 
specifying the absolute value of the angle between them. Hence we can see already that all n 
vectors must corotate. This can be confirmed further by the charges Rij calculated in (|3^) . 
They can be seen as two-dimensional determinants, fixing the area of the parallelogram defined 
by the above-mentioned two vectors. It then follows that the sign of the angle between these 
is determined, too. 

The other conserved quantities that involve only terms of first or lower order are the super- 



charges Qi and Qj, given in (pt). They couple the fermionic vectors with the bosonic vector: 
The charges Qi fix the scalar products, the constants Qi the determinants between {Xi, Pi) and 
{xq, U), so that all angles between these vectors, including their signs, are determined and con- 
stant in time. This means that once the bosonic motion is calculated and the fermionic initial 
data is given, we can read off the time evolution of all fermionic functions since their vectors 
are rigidly corotating with the bosonic one. 

Having now dealt with the lowest order bosonic and fermionic functions we have to ask how 
to interpret the higher order functions. Unfortunately, the geometric interpretation of these 
can not be seen as clearly. Instead we will interpret the higher order functions as variations of 
those of lower order. 

To start with, Xij as given in (^) has two terms coming from the two homogeneous solu- 
tions of the Xjj-equation with coefficients Cij and Eij and the term with coefficient Rij which 
constitutes a particular solution. We will now show that the first two terms can be written as 
variations of the lowest order motion xq with respect to the two free parameters t^ and Eq of 
that motion. 

The variation of Xo{t) with to is clearly proportional to the velocity Xq: 

Sxo 
dto 
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Figure 1: Bosonic (continuous) and fermionic (dashed) vectors, rigidly corotating 



which gives us the term with coefficient Cij. Next we vary the equation Hq = Eq: 

2xo5xo + 2UU'Sxo = 2 5Eq. 
Using equation (|26|) , dividing by Xq ^"^^ integrating between Iq and t we obtain 



to •'to 



hence 



and finally 



5xo 

Xo 

Sxq 

5En 



to 



/■* 1 
6Eo / —dt' 

J to Xq 



Xq / ) 

Jto Xq 



which is exactly the term with coefficient Eij in (|38D . So we end up with: 



6xn 



Xi 



—c 



5xn 



+ 



6to "^SEq ' 2Eo 



U. 



(61) 



In the same way we can analyse the third order fermion solutions P5| ) and Pq). Both have 
parts proportional to xq and U. But in addition there are two further terms which have to be 
interpreted as variations of these parts with to and Eq. The to- variations can be written down 
immediately: 

Sxq .. SU 

— OC Xo, — (XU Xq 

oiQ dto 

and yield the terms with coefficients C[ijlk] and C[ijrk]. The E'q- variations give: 



6xn 



d 6xo 



5Eo dt 5Eo 



^0 / ^dt' + 1 1^ = U'^ = U'xo f ^dt' 

JtQ Xq Xq OH/ci OH/f) Jtn Xn 



Xq' 5Eo 



i 1 

to Xg 
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They explain the first part of the terms with coefficients E^ijlk] and E^^rky, note that the second 
part can be absorbed into the Ujk and rjjfc-terms. 

So the vectors formed by Xijk and pijk consist of three parts: A vector rigidly fixed to the 
familiar first order fermionic vectors and two parts which are proportional to the variation of 
these vectors with and Eq. 

Finally we want to explain the parts of the highest order bosonic function Xi234- The 
first three terms of ( |1D| ) are familiar: They are the variations of Xq with to and Eq and the 
inhomogeneity term proportional to U. The next three terms can be explained as second 
variations. So xq gives the second variation of xq with respect to to] for the term with coefficient 
C[i2-£'34] we have already shown that it is the time derivative of hence it is proportional to 
j^^. So we only have to deal with the £'[12 £"34] -component. We find that 

= 7^ ^0 / ^dt = -— / —dt + XoYTT / —dt ■ 
SE^ SEo \ Jto J 5Eq Jto xg 5Eq Jto xg 

The last term can be rewritten as 

/■* 1 / 1 /■* 1 

= -2xo / —dt' + xo — / —dt' 

Jto Xq J^O Xq 

Combining this with the expanded first term 

Xo ( / ^dt'] + ^ / -^dt' 

\Jto Xq j Xq Jto Xq 

we find exactly the component in (^0]) with coefficient £'[12-^34]. 

The last two terms in (^0|) are easy again: They are given by the two first variations of the 
inhomogeneity term with to and i?o, which we have already calculated above. Note that 
the term — with coefficient i?[i2-E'34] can be absorbed into the -R1234 term. 

In summary: 

1. The lowest order bosonic function xq describes the one-dimensional motion of a point 
particle with energy Eq in the potential |?7^(xo). It is completely unaffected by all higher 
order bosonic and all fermionic functions and not influenced by any higher order charges. 
There are two free parameters involved, an initial time to and the energy Eq. 

2. The first order fermionic functions comprise two-dimensional vectors (Aj, pi). The lengths 
of these vectors and the angles between them are fixed by charges arising from symmetry 
under change of basis in the Grassmann algebra and from the second order charge 
and so are constant in time, i.e. all fermionic vectors rigidly corotate. Their motion is 
in turn determined by the bosonic vector (xq, ^) which has constant squared length 2£'o 
and is rigidly coupled to the fermion vectors by the supercharges Qi and Qi. For every 
generator there is a pair of free parameters which specify the length and the phase of the 
corresponding vector. 
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3. The second order bosonic function xtj has terms proportional to the (first) variation of 
Xo with to and Eq. There is one additional term proportional to U stemming from the 
inhomogeneity of the Xjj-equation of motion and ultimately from the Yukawa interaction 
term. 

4. The third order fermionic functions again comprise two-dimensional vectors {Xijk, Pijk)- 
They can be divided into three additive parts. The first part is rigidly rotating with the 
first order fermionic vectors, the second and third parts are proportional to the variations 
of this first order vectorial motion with initial time to and Eq. 

5. Finally, the fourth order bosonic function involves all first and second order variations 
of Xo with respect to to and Eq plus a further term proportional to U and its first order 
variations. 



9 Oscillatory Motion 

In this section we apply our results of the previous section to oscillatory motion, i.e. we assume 
that the lowest order bosonic function xo is periodic with period T. We have found already 
that in this case the first order fermionic motion is described by two-dimensional vectors that 
either continuously wind around a circle or oscillate between two symmetrically lying points on 
it, depending on the number of sign changes of U. 

When we look at the next bosonic level, i.e. the functions Xij, we can see immediately that 
the first part of the solution (pSj), which is proportional to xq, is a periodic function with the 
same period T. The last term proportional to U{xQ{t)) is periodic, too, but it can have period 

if the potential function U{xq) has reflection symmetry, i.e. if U{a + x) = U{a — x) for some 

constant a. Let therefore T, which is either T or |T, denote the period of [/(xo(t)). The most 
interesting term is the remaining second one containing the integral. We will analyse this term 
from a mathematical point of view using Floquet theory but before we do that we would like 
to understand it from a more physical perspective. 

From equation (^) in the previous section we know that the integral term can be interpreted 
as the variation of the solution xq with respect to the energy Eq. Because the variation is 
infinitesimally small, we can make the assumption that the functional form of the motion 
remains unchanged. We treat the motion as determined by only two parameters, the period T 
and a characteristic amplitude A, defined e.g. as the distance between the two turning points 
of the motion. In the generic case T and A will be related by a non-trivial function T{A), 
therefore a change in energy means not only a change in the amplitude A but also in the period 
T, so: 

dxQ dxQ dA ( dxQ ^ Oxq dT\ dA 



dEo dA dEo \ dA dT dA J dEo 

^] ~ T- 



Since xq is a T-periodic function we can use its Fourier series CjC*'^-'*, where Uj = %j and 



the coefficients Cj are regarded as functions of A, to find: 

BT T ' J J 

3 

Thus 

dxn , , / dxn , , Xn(t) dT\ dA 



d^o V T dA dEo 
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Both ^ and Xq are T-periodic functions, but itself is clearly not since the second term 
is /mear-periodic, i.e. the product of the linear term t with a periodic function. The extra 
linear-periodic term diverges with time, thus seemingly spoiling the oscillatory nature of the 
solution. However, this problem only arises because the period of the oscillation depends on the 
amplitude A and thus ultimately on the energy Eq. Expressing a solution with slightly changed 
period T + dT in terms of T-periodic functions inevitably leads to non-periodic terms. The 
effect is well-known from celestial mechanics where linear-periodic functions appear as secular 
terms in the study of stability problems of planetary orbits (see e.g. 0). 
There is one case where no linear-periodic term occurs, and this is when 

dT 

i.e. when the period is independent of the amplitude. This is true for the harmonic oscillator, 
and consequently there is no non-periodic term in equation (^). 

The existence of non-periodic terms in the bosonic function be derived in a math- 

ematically more stringent way from Floquet theory, the theory of linear differential equations 
with periodic coefficients. Equation (p7D, with Xipj — piXj replaced by the constant Rij ac- 
cording to (p4D, determines the motion of Xij. Because we know that the particular solution 
proportional to U is periodic it suffices to treat the homogeneous equation which can be written 
in the form 

^ij +p{t)xij = 0, (62) 
where p is a T-periodic function. This allows us to apply Floquet theory which states that [|]: 

1. There exists a non-zero constant a, called the characteristic multiplier, and a non-trivial 
solution Xij{t) such that 

Xij{t + T) = axij{t), (63) 

from which one deduces 



2. There exist linearly independent solutions Xij^i and Xij^2 to (|6^), such that either 



X,,- i(t) = e'"i*Pi(t), x,,,2it) = e™^*P2(t) (64) 

or 

X,,- i(t) = e™*Pi(t), x^,,2{t) = e'"*(tPi(t) + P2(t)), (65) 

where in both cases Pi, P2 are T-periodic functions and mi, 1712, m, called characteristic 
exponents, are - not necessarily distinct - constants. 

Whether the solutions take the form ( |5^ or ( p3| ) depends on whether there are two independent 
solutions of (|62D with the property (^) or just one. If we denote by ^1,^2 the two linearly 
independent solutions of (|62D with 

Xi{0) = l,Xi(0) = 0,^2(0) = 0,X2(0) = 1, (to = for simplicity) 

and by M{t) the matrix 

/ X^it) X2{t) \ 
I X,{t) Mt) j ' 
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we get a solution of type (|64D if M{T) is diagonalizable and a solution of type (^) if M{T) 
has Jordan normal form. 

To determine the characteristic exponents nii it is useful to notice that our equation is not 
an arbitrary differential equation with periodic coefficients but an example of Hill's equation 
which takes the form 

F{t)xij + F'{t)xij + G{t)xij = 0, 

F and G being T-periodic functions. For equations of Hill's type the Floquet multipliers are 
completely determined by the trace D of the fundamental matrix M{T). This trace can be 
determined indirectly, using the fact that we know one periodic solution of (|6^) , namely Xq, 
which has either period T or 2T, depending on whether U{xo) is reflection-symmetric or not. 
Now we can use another theorem of Floquet theory which states: 

1. Hill's equation has non-trivial solutions with period T if and only if D = 2. Then either 
TTT-i = m2 = m = and 

or m = and 

XijA^) = Xij,2it) = tPi{t) + P^it). 

2. Hill's equation has non-trivial solutions with period 2T if and only \i D = —2. Then 
either mi = m2 = ^ and 

ITT ^ iiT ^ 

or m = ^ and 

x^J,l{t) = et*Pi(t), Xij,2{t) = ef\tP^{t) + P2{t)). 

Using these theorems we find that the second independent solution ( ]61| ) has to be either pe- 
riodic or linear-periodic. The period is in both cases T, the period of the first independent 
solution xq. This means that we can formally confirm our earlier result derived from physical 
arguments. Additionally, we have that any solution must be bounded in finite intervals. This is 
not immediately obvious from the explicit form of (|6l|), which could be singular when xq = 0. 

Before we proceed to analyse a particular potential function which admits oscillatory motion 
we want to comment briefly on the higher order bosonic and fermionic quantities Xijk,Pijk and 
3^1234- From equations P5| ) and (^6[) we can see that all terms that appear in the third order 
fermionic solution are either periodic or linear-periodic. For the term 

/■* 1 1 
^0 / —dt' + — 

Jto Xq Xq 

this follows e.g. from the fact that the derivative of a linear-periodic term is again linear- 
periodic. We can expect this behaviour physically since Xijk and pijk are given by the first 
variation of the motion of the first order fermionic vectors with energy and initial time. (There 
are further periodic contributions to Xijk and pijk with coefficients kjk and Vij^.) 

Most terms in the fourth order bosonic solution X1234 can be easily seen to be periodic 
or linear-periodic in the same manner. However, there are the three terms which have been 
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interpreted as the second variation of Xq with energy Eq. These should contain a quadratic 
term in t, because 



dEl 



d'^xo d'^xo dT d'^XQ 




+ 



dxp d^T 
dfdA^ 



'dAT 

dEa 



+ 



dxo dxo dT 
'dA ^ WdA 



dEl 



and 



d'^xo 



d_ 
df 



Xot 



1 



2xot + Xot' 



all other terms being (linear-)periodic. Again, the quadratic term does not occur when period 
T and amplitude A are independent, i.e. when ^ = 0. This is confirmed by the harmonic 
oscillator example where X1234 is completely periodic. 

We can also understand the existence of a quadratic-periodic term directly from the 
solution ( ^OD which contains the term 

Xo( f -^dt' 

\Jto Xq 

We know already that the integral will yield a term proportional to t; this is squared to give 
a quadratic term in t, and then multiplied by the periodic function xq. This quadratic term 
cannot cancel out because the two remaining terms with coefficient i?[i2-E'34] iii (0) most 
linear-periodic. 



9.1 The hyperbolic potential U{x) = ctanhkx 

We illustrate our results by the hyperbolic potential U{x) = ctanhfcx, already discussed in 
section |4.2| for general Grassmann algebra. Since we are only interested in oscillatory motion 



here, we make the assumption that > 2Eq. 

Inverting the result obtained from (BB) we can calculate that 



k 

with 



Xo{t) = i^arcsinh ( ^^5^ sin(a;(t-to) + , (66) 



u = k\/ (? — 2Eq , K = arcsin (^y — 1 sinh kxaito] 
So Xq is indeed a periodic function. For xq and U we immediately find that 



^o{t) = ± (67) 

, , ctany 
U{t) = ± / (68) 

where we have abbreviated {uj{t—to) + K,) to y. Therefore the motion of the bosonic vector ( ^9]) is 
a non-uniform rotation, which confirms our comments in section ^ since the potential function 
U changes sign only once. All fermionic vectors corotate with the bosonic one. 
To determine Xij we need the integral 



—dt' = — . -uj{t-to) + — 

*o x'o kJ{c^-2Eoy \ 2E, 
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we then obtain for the i^jj-related part of Xj,-: 



± 



Ej, 



c 

2Eo 



tany 



2Eo 



sec^ y — 1 



^sec^y-l^ 



(69) 



This result clearly shows a linear periodic term which means that the period-amplitude relation 
is non-trivial. In fact, period and amplitude of Xo{t), the latter defined as distance between the 
two turning points, can be easily calculated using (p6|): 



T 



27T 



, A = — arcsinh^ 



2Eo 



2En 



Thus 



T{A) 



27r , kA 
— — - cosh — 
k\c\ 2 



and ^ 7^ 0. We point out that all terms in the solution are clearly bounded, which is not the 
case for the integral itself. 

Figure ^ gives an example of both the bosonic quantities Xq and Xij and the fermionic 
quantities Aj and pi. 




Figure 2: xq (continuous), Xij (dotted), Aj (long dashes) and pi (short dashes) as functions of time 
for c = k = I, xq{0) = Cij =0, Eq = j, Eij = i and Qi = Qi = ^ 



We do not present the higher order solutions in detail here since their explicit form, although 
not complicated to find, is lengthy and does not provide new insights. The only interesting 
exception is the part of X1234 which has been characterized as a second variation with energy. 
This part can be calculated to be 



—EuoE' 



[12-C/34] 



+ (t-to) 



(a'^ — Sa^ + Ga) sinysec^ y + (2a^— 4a) siny secy 



A;(c2-2Eo)y(asec2y-l)3 
;-2a2 + a) sec^ y+l 

+ ~ 'foj 



ka 



sm y sec y 



c2-2Eo)^/(asec2y-l) 



c^-2Eq ^/(asec^y-l) 



where we have abbreviated to a. Again, there are two noteworthy points here: Firstly, the 
occurence of a quadratic-periodic term, which has been predicted in section ^ and secondly the 
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complete boundedness of the whole solution, which is not immediately obvious from the func- 
tional form (HDf), especially regarding the term where the velocity appears in the denominator. 



Next we want to show that the two solution methods presented in this paper deliver the 
same result. Therefore we have to decompose our previous solution, the first of equations (|2T|) , 
from section [4.2| . Recall that any real analytical function / is extended to a Grassmann- valued 
function by the formula 

OO 1 



n=0 



where Zb and Zg are body and soul of the Grassmann number z, respectively 0. To keep the 
calculation short we will expand here only up to second order, i.e. the following equations have 
to be understood modulo a component proportional to ^1234- 
We begin with the frequency u: 

kE, 



u = kJ{c^-2Eo)-2E,,^,, = k^c^-2E, - 

^ \I& — IEq 

Then 



sin(a;(t — to) + = sin(A;\/c^ — 2£'o(t — to) + ^o 
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Since 



l2E+2ckR I 2Eo I c^Ea + ckRiJc^ - 2E^ 
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'2Eo{c^-2Eor 



we find 



l2E+2ckR , , , , 2Eo 2Eo 

sm{u{t-to) + K) = d smy + 



c2-2E ' ' ' V c^-2Eo ^ V c2-2^, 



-^ck sin y + Kij cos y 
2Ea y V y 



+Ei. I sin y - cos y 







where again y = k\J (? — 2Eo{t — to) + kq. This result equals sinhfcx; decomposition yields: 

sinhfcx = sinh(A;a;o + kxij^ij) = sinhA;xo + kxij{cosh. kxQ)C,ij , 
thus by comparison 

sinh kxo = J — — sin y, 



c^-2E< 







which gives the same solution for xq as ([66|). Similarly, we can read off the result for kxij cosh kx^ 
calcu 

Xij it) 



and then calculate xij 



R'l j ctan?/ Ky^j 



^^°\/&sec2|/-l k y^^sec2y-l 

2^ tany 1 t-to 



+E.. 



k{c^-2Eo) ^^sec^y-l Vc^-2^o^^sec2y-l^ 
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The term with coefficient Rij can be readily identified as the potential function U{t), derived 
in (p^). The second term with coefficient Kij equals the velocity xq, which we have calculated 
in ([67|), when we make the identification = —uocij. Finally, the term proportional to Eij is 
identical with our result (^Uj). 

So the solution found by using our layer-by- layer approach can also be obtained by the de- 
composition of the full Grassmann solution calculated through direct integration in Grassmann 
space. We have demonstrated this feature here using the bosonic quantity x{t) but we could 
equally well have chosen one of the fermionic quantities ip+{t) or i^-{t). This means that our 
two solution methods are compatible with each other. 

10 The Zero Energy Solutions 

We have already mentioned in sections ^ and ^ that the case where the energy Eq, the body of 
the full Grassmann energy E, equals zero has to be treated differently. All the solutions to the 
equations of motion (p6D-(|30D were given under the restriction i^o 7^ so that we cannot use 
our previous results. It will turn out, however, that for a finitely generated Grassmann algebra 
we can derive all solutions in explicit functional form for arbitrary potential. 

Starting with equation (|3TD we can see that if i?o = both xq and U (xq) have to vanish: 

xo = 0, U (xq) = 0. 

This means that the particle stays permanently at rest at a minimum xo^mm of the potential U'^. 
One could assume at this point that all higher order bosonic and fermionic functions are trivial 
as well, but as we shall soon see this is far from being necessary, in contradiction to 0. 

Because xo(t) is constant, all the (spatial) derivatives of U are constant functions of time 
too; we denote their values by 

U'{xq) = ki, U"{xo) = /C2, U"'{xo) = ^3, . . . ; 

we assume in the following that fci 7^ 0. 

When we now look at the first order fermionic equations (p9D for Aj and pi, we find that 
they are easily solved by the harmonic oscillator ansatz: 

A, = li cos ki{t — to) — ri sin ki{t — to) (70) 
Pi = k sin ki{t - to) + Ti cos ki{t - to). (71) 

Unlike the true harmonic oscillator, the constants k and are, however, not linked to the first 
order supercharges Qi and Qi which vanish completely as can be verified from (|4TD. 

To understand the first order fermionic motion we assume that the energy Eq is small but 
non-zero, restricting the particle motion to a small neighbourhood of the stability point Xo,mm 
where the potential function U is approximately linear. The result is an almost harmonic 
oscillation with frequency ki, mirrored by the fermionic quantities - as follows from (|48| ) and 
(^). When we now let Eq approach zero, this fermionic motion seems to diverge as one can 
see from the formulae for the coefficients Li and Ri in (^Dj) and (0). However, there is a subtle 
point here: While the supercharges Qi and Qi are arbitrarily chosen constants for £'0 7^ 0, they 
have to vanish for Eq = 0. Thus to avoid any discontinuities they have to smoothly approach 
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zero as the energy decreases. Exactly how they approach zero finally determines which value 
the harmonic oscillator coefficients Li and Ri take in the limit Eq = 0, which justifies the two 
remaining degrees of freedom in our solution, li and rj. 

Next we want to analyse the second order bosonic quantity Xij. Because U' and U" are 
constants the equation of motion ( pTl) describes formally a harmonic oscillator with frequency 
ki subject to a constant external force k2Rij, so the solution (for ki ^ 0) is 

^2 

Xij = Vij cos ki{t - to) + Vij sin ki{t - to) + ja^iji i^'^) 

ki 

where Vij and Vij are integration constants. To interpret this result we recall that in case Eq 
the homogeneous part of Xij consists of two terms, the two variational derivatives of xq with 
respect to Eo and to- When Eq is zero, a small change in energy will result in oscillatory 
motion around the stability point a;o,mm- This explains the functional form and one of the free 
parameters of the solution (|7^). The second parameter, though, cannot be connected to the 
variation with to anymore since this variation is zero for a constant Xq. 

An interesting observation is connected to the second order energy which can be calculated 
from (|32|) to be Eij = —kiRij. Thus the second order energy is not independent any more but 
determined by the four first order constants li, rj, Ij, rj. Note the fact that Eij is not connected 
to the parameters Vij and Vij. 



The third order fermion equations ( PU] ) are again equations describing a harmonic oscillator 
with a driving term. The homogeneous part of the solution looks therefore the same as ([70|) 
and (|7TD with integration constants kjk and Vijk. To find the particular solution we have to 
investigate the driving terms —k2X\ijPk] and k2X[ij\k] further. Using equations (^), ([7T|) and 
(0) we find 

X[ijPk] = Cijk + Di^ijk sin 2ki{t - to) - D2,ijk cos 2ki{t - to) 
X[ij\k] = Cijk + D2,ijk sin 2ki(t- to) + Di^ijk cos 2ki(t- to), 

where Cijk, Cijk, Di ijk and -D2,ijfc are constants built from li, ri, Vij and Vij. So the driving terms 
are not constant but have an oscillating part that oscillates with double the basic frequency ki. 
The particular solutions are thus 

Kjk,part. = —j^Cijk + -^Di ijk cos2/ci(t — to) + -^D2^ijk sin2A;i(t — to) 

Pijk,part. = ~Y'Cijk + 7— -Di^jjfc sin 2/ci(t — to) — -^D2^ijkCos2ki{t — to). 

The fermionic vectors {Xijk, Pijk) consist of three parts: A uniform rotation with natural fre- 
quency ki, a forced uniform rotation with double this frequency that we can interpret as the 
first excited mode and a constant shift away from the origin. Note that there is no double 
frequency term if U"{xo) = 0, and in particular for the harmonic oscillator potential. This 
agrees with our earlier result in section |6.1| . 



Calculating the supercharges Qijk and Qijk from (|l|) we find that Qijk = 2kiCijk and 
Qijk = —2kiCijk which gives us a neat interpretation of the displacement constants Cijk and 
Cijk- It further means that the third order supercharges are not related to the free parameters 
lijk and rijk as is the case for Eo but instead are completely determined by first and 
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second order parameters. This is the same phenomenon aheady encountered for the second 
order energy Eij and again it means that we should not think of Qijk and Qijk as independent 
constants when Eq = 0. 

Finally we have to investigate the fourth order bosonic quantity xusa. The equation of 
motion ( ]28|) simplifies in the case i?o = to 

3 

3^1234 + ^13^1234 = --klk2X[12Xu] + A;2-Rl234 + ^3a^[12^34], 

SO it describes a harmonic oscillator with three driving terms. The homogeneous part is analo- 
gous to ( [7^ ) with two integration constants C1234 and C1234; the comments made about Xij apply 
here, too. Analyzing the three driving terms we first note that ^2-^1234 specifies a constant 
external force, playing the same role as k2Rij for the variables Xij. The first and third terms 
can be calculated using the explicit formulae for Xij, resulting in 

3 3 
--kik2X[i2Xu] + k3X[uR34] = -^kik2 (Ci234 + 1)1234 COS 2/ci (t - to) + Di234 siu 2A;i(t-to 

k^\ / 

-3^1 (Fi234 cosfci(f-to) + Fi234 sin/i;i(t-to 

where C1234, -D1234, -D1234, -F1234 and F1234 are fixed constants built from Vij, Vij and Rij. So 
besides a further constant term we find four oscillating ones, two with frequency ki and two 
with frequency 2ki. The latter ones will cause forced harmonic motion with double the natural 
frequency as we found for the third order fermionic terms. New, however, are the oscillating 
terms which have the same frequency as the homogenous solution. This means that the os- 
cillator is driven in resonance. Therefore the solution is not bounded anymore but includes 
linear-periodic terms like tsinki{t—to) and tcoski{t—to). Combining all contributions we find 
the particular solution 

Xl234,part. = y (^^J^-jCuSi) + [D 122,4 COS 2ki{t - to) + I)i234 siu 2A;i(t-to; 





^2 k \ 

3^--^) (Fi234 icosA;i(t-to) - F1234 ^ sin /i;i(t-to) 



Thus the fourth order bosonic solution consists of harmonic oscillation with the natural fre- 
quency ki, an excited mode with twice this frequency, linear-periodic motion and a constant 
shift. 

Especially interesting are the linear-periodic terms. They can be explained physically as 
in section |^, although we have to see them as the second variation with energy here: When 
the particle stays at rest in a potential minimum the first variation gives rise to harmonic 
oscillation, independently of the shape of the potential function U (as long as /ci 7^ 0). Only 
when we look at the second order variation do non-harmonic terms come into play: The period 
of oscillation will in general depend on the amplitude and the energy, therefore a variation in 
energy results in a change of period; expressing the new solution in terms of the old period 
then generally leads to secular, non-periodic terms. 

Notice the fascinating fact that although the physical interpretation is similar, the way the 
linear-periodic terms arise mathematically is completely different: For non-zero energy oscilla- 
tions they come as solutions of a homogenous differential equation with periodic coefficients; 
in the zero energy case they are particular solutions to an inhomogenous differential equation 
with constant coefficients. 
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11 Discussion 



In this work we have analysed a supersymmetric mechanical model from two different view- 
points: Either we make no specification with regard to the nature of the underlying Grassmann 
algebra B. Then for a range of potential functions U the model can be explicitly solved. Or we 
regard B as finitely generated thus reducing all quantities to a set of real functions and their 
interrelationships. Then we are able to solve the system completely, without any restrictions 
on U. The methods have been shown to be compatible with each other here. 

An open question with our first point of view is what meaning we can give to the full, i.e. 
non-decomposed Grassmann solutions x{t),ip+(t) and il^-{t). Interpretation may be aided by 
our second approach, where we have shown that the component solutions include the purely 
classical motion dependent on Eq and all its variations with respect to this (real) constant. 
Apparently, supersymmetric dynamics captures information about a whole range of energies of 
the mechanical system, so in a sense we can say that the Grassmann energy E corresponds to 
some fuzzy classical energy. This suggests possibly that supersymmetric classical dynamics is 
closely related to the quantum dynamics. A study of the quantized version of our model would 
therefore be an important step to complement this research. 

We also hope to apply our methods to the full, i.e. space- and time-dependent field theory. 
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